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"Method for decoding linear space -time codes in a multiple -antenna 
wireless transmission system and decoder for implementing such a 

method. " 



The present invention relates to a method for decoding linear 
space-time codes in a multiple-antenna wireless transmission system. 
It also relates to a decoder implementing this method. 

The invention is found to be particularly useful in the field 
of radio transmission or broadcasting of digital data, or sampled 
analog data, in particular in the case of transmission with mobiles 
or more generally, in the case of a local or non-local wireless 
network. More precisely, the invention can be used in particular when 
high-throughput wireless transmissions are to be made. A first 
category of use relates to cellular communication systems with 
mobiles such as the UMTS for example. A second category of use 
relates to local wireless networks. A third category of use is that 
of future ad hoc networks. 

Existing and future wireless communication systems require a 
high quality of transmission for constantly increasing throughputs. 
In fact, these systems have as an objective the establishing of 
services including high quality data and video transmission channels. 

There are numerous constraints linked with multi-path 
propagation. This multiple propagation is due to the fact that a 
radio signal reaches an antenna by several paths via reflections on 
the ionosphere and on terrestrial objects such as mountains and 
buildings for example. The effects of this multiple propagation are 
constructive and destructive interferences as well as phase shifts of 
the signals transmitted. 

In order to develop a reliable and robust communication system, 
space, time and frequency diversity techniques are used. Space 
diversity consists in particular of an arrangement of several 
transmit and/or receive antennae. When time diversity is added, a 
space-time diversity is created, 



requiring space-time coding on transmission, as well as space-time 
decoding on reception. It is known that space-time codes make it 
possible to increase the throughputs of information transmitted for 
an equal quality of service. 

High-efficiency linear space-time codes have already been 
decoded by decoding algorithms which can be classified in two 
families : 

interference cancellation algorithms such as BLAST ("Bell 
Laboratories Layered Space-Time") which have poor 
performances in terms of error rates as a function of the 
signal-to-noise ratio. 

maximum likelihood algorithms which are optimal in terms of 
performance, but with a high installation complexity such 
as the sphere decoder. 
The document US 6 178 196 (SESHADRI et al . ) is known, which 
describes a multiple antenna system combining interference 
cancellation decoding with maximum likelihood decoding. The system 
described has the purpose of isolating the signals originating from a 
given station, assuming that the symbols of the signals transmitted 
by another station are correctly determined during a preliminary 
estimation . 

An optimal decoder for high-yield linear space-time codes has 
been found by Damen, Chkeif and Belfiore (O. Damen, A. Chkeif and J.- 
C. Belfiore, "Lattice Code Decoder for Space-Time Codes," IEEE 
Communications Letters, May 2000) . Its performances exceed those 
possible with the "BLAST" decoder (G.D. Golden, G.J. Foschini, R.A. 
Valenzuela, P.W. Wolniasky, "Detection Algorithm and Initial 
Laboratory Results using the V-BLAST Space-Time Communication 
Architecture," IEEE Electronics Letters, January 7th 1999). However, 
it suffers from three problems closely linked to its implementation: 

1. It is of a fairly high complexity, which makes its high throughput 
implementation very problematical. 

2. Its complexity strongly depends on the radius of the sphere 
chosen. The code word decoded must in fact be found in this sphere if 
the decoding is not to fail, and moreover, there must not be many 
other code words to be found in this same sphere, 
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as otherwise the decoding complexity can further increase very 
considerably . 

3. Finally, the decoding complexity depends very strongly on the 
signal-to-noise ratio and also the right choice of sphere. The sphere 
choice algorithms are themselves fairly complex. It is therefore seen 
to be necessary to find another decoding algorithm which does not 
suffer from this kind of problems. 

The present invention has the purpose of developing a new 
linear space-time decoder realizing a compromise between the 
interference cancellation technique and the maximum likelihood 
technique. The invention aims at performances very close to the 
maximum likelihood technique with simplified implementation compared 
with that of the interference cancellation algorithms. 

The abovement ioned objectives are achieved with an iterative 
method for decoding a set of N sampled signals in a space-time 
communication system with M transmit antennae and N receive antennae. 
The N signals are processed by intervals of time T corresponding to 
the time length of the linear space-time code associated with the 
transmission signals. In each time interval T, the sampled signals 
received are represented in the form of a signal matrix Y. 



can also be expressed in the form Y = H X + W (1) 

in which X is the space-time code word transmitted and is represented 
by a rectangular matrix with M lines and T columns, H is the channel 
matrix and is represented by a rectangular matrix with N lines and M 
columns, and finally W is noise, represented by a rectangular matrix 
with N lines and T columns. 

It is therefore seen that the line changes correspond to the 
antenna changes whereas the column changes correspond to changes in 
the sampling times. In this model, the coeffecients of the matrix X 
depend linearly on the information symbols to be transmitted, i.e. 
the coefficients of the matrix transmitted, Xij with i ranging from 1 
to M and j ranging from 1 to T are written in the form 
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with s pq being the information symbols and gij pq being coefficients 
which depend on the code chosen. Thus, any linear space-time code, 
i.e. such that the words transmitted X have their coefficients x ijf 
which can be expressed in the preceding form, can be decoded by the 
method according to the invention. 

A description will now be given of a vectorized model on which 
the method according to the invention can be used. Instead of working 
with matrices, it is possible to vectorize the expression of the 
signal received in equation (1) and write 



15 



Y = vec (Y) = H X + W (2) 

With x = vec (X) = G S, G has the dimension (MT,MT) 



And 



H = 



H 



0 



H 



20 H is a matrix (NT, MT ) having on the diagonal T times the matrix 

H, the other coefficients being zero. 

Equation (2) is then equivalent to equation (1) and constitutes 
its vectorized version. It is on this version that the decoder of the 
invention will be used. It is now sufficient to define the extended 

25 channel matrix 

C = H G 



which will be used hereafter. The vectorized version of the received 
30 signal can therefore be rewritten 



Y = C S + W 
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where Y is a column vector with NT components, C is a matrix NTxMT, S 
is the symbols column vector with MT components and W is the noise 
column vector with NT components. 

The iterative decoding method according to the present 
5 invention makes it possible to obtain an estimation of the symbols of 
the signals transmitted. This method is used in particular for N 
greater than or equal to M. According to the invention, it comprises 
the following steps: 

Pre-processing of the vector Y in order to maximize the 

10 signal to noise+interf erence ratio in order to obtain a signal r 1 . 

- Subtraction from the signal 7' of a signal Z l by means of a 

subtractor, the signal Z* being obtained by reconstruction post- 
processing of the interference between symbols from the symbols 
estimated during the preceding iteration; 

15 - Detection of the signal generated by the subtractor in order 

to obtain, for the iteration in progress, an estimation of the 
symbols of the signals transmitted. 

The pre-processing step can be carried out by operating a 
matrix multiplication between the signal vector Y and a matrix B, the 

20 matrix B being updated at each iteration. 

The post-processing step can also be carried out by operating a 
matrix multiplication between the vector of the symbols estimated 
during the preceding iteration and a matrix D, the matrix D being 
updated at each iteration. 

25 According to the invention, for each iteration, a standardized 

correlation coefficient p is calculated, the updating of a matrix 
being achieved by determining new coefficients of the matrix as a 
function of the correlation coefficient obtained for the preceding 
iteration . 

30 According to a preferred embodiment of the invention, the N 

signals are processed by time intervals T corresponding to the time 
length of the linear space-time code associated with the transmission 
signals, and the pre-processing step involves the matrix B in order 
to maximize the signal to noise+ inter f erence ratio, the transfer 

35 function of which is: 
/ \ 



B = Di °s 



i 



E s ) 



C H V f 

\<i<MT 
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Fith V F = 



l-A 2 -i 



No 



CC H +Id N 



; j( =diag(c» V< C); 
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€ : iteration index; p: standardized correlation coefficient between 
the real symbols and the estimated symbols; N 0 : noise variance; Es : 
mean energy of a symbol; C: extended channel matrix. 

Similarly, the post-processing step can involve a matrix D for 
the reconstruction of the interference between symbols, the transfer 
function of which is: 



E>=B<.Cp^-Diag 



1 



\<i<MT 



Before the first iteration, we have no information on the 
symbols. The matrix B has the role of maximizing the signal/ (noise + 
15 interference) ratio. The vector z is void. A representation of the 
vector of the symbols S can be such that the k th symbol is: 



S k (received) = S k (transmitted) + / M Q. : S '.- (transmitted) + noise 

20 At the second iteration, the matrix B still maximizes the SINR 

ratio. The matrix D will mimic the interference between symbols, i.e. 

^1 CLi S: (transmitted) during the preceding iteration at the level of the 

signal r leaving the matrix B. The subtractor makes it possible to 
subtract this interference. 
25 At the last iteration, it is assumed that the symbols are 

correctly estimated, i.e. that D makes it possible to reconstruct all 
the interference, such that it is estimated that: 
S k (received) = S k ( transmitted) + noise 

The invention is a decoder which can be adapted to linear 
30 space-time codes, whatever they are, i.e. it makes it possible to 
decode any code such that the sequences transmitted are written as a 
linear combination of the information symbols. 
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It can be noted that the matrices B and D depend on the 
correlation coefficient which is different for each iteration. 
According to an advantageous characteristic of the invention, in 
order to determine the correlation coefficient p e , at each iteration: 
5 - the signal to interference ratio SINR is calculated at the 

threshold detector's input using the following formula: 



SINR' = 



-1 
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and — 
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- the symbol error probability Pr, for example at the threshold 
detector's input, is calculated from the signal to interference ratio 
SINR'; and 

20 the correlation coefficient p' is then calculated from the 

symbol error probability Pr . In order to do this, it is possible to 
use a standard formula producing the correlation coefficient as a 
function of the symbol error probability, this formula depending on 
the modulation used in the transmission. 

25 This results in the most precise possible approximation to the 

correlation coefficient p e . 

In subsequent experiments, four iterations sufficed to obtain 
very good results. But it is also possible to define a minimum value 
(threshold value) of the correlation coefficient for which the 

30 iterations are interrupted. This coefficient is essentially a 
function of H, which can be estimated in standard manner, and of the 
variance N 0 . Consequently, all the values of the correlation 
coefficient p, as well as the values of the matrices B and D (for all 
of the iterations) can be calculated before the first iteration. They 

35 can therefore be stored 
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before the first iteration in a memory in order to then be brought 
out for each iteration. 

Preferably, it is assumed that p° = 0 . 

Moreover, in order to calculate the symbol error probability 
5 Pr, it can be assumed that the total noise is Gaussian and that it is 
possible to use the formula corresponding to the constellation of a 



Shift Keying") system or to use tables indicating the error 

probability as a function of the signal to noise ratio. In fact, 
10 depending on the complexity of the symbol error probability it can be 

useful to directly tabulate the formula. 

Advantageously, in order to calculate the correlation 

coefficient p f from the symbol error probability Pr, it can be 

assumed that when there is an error, the threshold detector detects 
15 one of the closest neighbours to the symbol transmitted. 

By way of example, at the final iteration, it is possible to 

introduce the signal leaving the subtractor into a soft-entry 

decoder . 



20 invention, the information symbols can be elements of a constellation 
originating from a quadrature amplitude modulation, or QAM. 

According to another feature of the invention, a space-time 
decoder is proposed for decoding a signal vector Y obtained from N 
sampled signals in a space-time communication system with M transmit 

25 antennae and N receive antennae, with N greater than or equal to M, 
with a view to obtaining an estimation of the symbols of the signals 
transmitted. According to the invention, this decoder comprises: 



linear modulation, for example 




the BPSK ("Binary Phase 



According to a preferred method of implementation of the 



30 



- a pre-processing module of the vector Y for maximizing the 
signal to noise+interf erence ratio in order to obtain a signal 7* , 

- a subtractor for subtracting a signal z 1 from the signal 7* , 
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a post-processing module for the reconstruction of the 
interference between symbols from the symbols estimated during the 

preceding iteration in order to generate the signal Z , 

- a threshold detector for detecting the signal generated by 
5 the subtractor in order to obtain, for the iteration in progress, an 

estimation of the symbols of the signals transmitted. 

These pre-processing and post-processing modules can be 
matrices, B and D, according to the formulae indicated previously. 

Other advantages and characteristics of the invention will 
10 become apparent on examining the detailed description of an 
implementation embodiment which is in no way limitative, and the 
attached drawings in which: 

Figure 1 is a diagram illustrating some elements of a 
transmission system within a transmitter and a receiver, the space- 
15 time decoder according to the invention being integrated into the 
receiver; 

- Figure 2 is a general diagram illustrating the architecture 
of the space-time decoder according to the invention; and 

Figure 3 is a general block diagram according to the 
20 invention. 

Figure 1 represents a transmitter 1 provided with a plurality 
of antennae 7. The transmitter 1 comprises in particular, upstream of 
the antennae, an error correcting encoder 3 followed by a linear 
space-time encoder 4. According to an advantageous characteristic of 

25 the invention, the decoder according to the invention is capable of 
being applied to any linear space-time code, i.e. codes such as the 
sequences transmitted are written as a linear combination of the 
information symbols. 

The signals transmitted by the antennae 7 are picked up by a 

30 plurality of antennae 8 within a receiver 2. The signals received 
undergo processing within a space-time decoder 5 according to the 
invention in order to estimate the information symbols of the signals 
transmitted. The space-time decoder 5 has a soft output to which a 
soft-input decoder is connected, for decoding error correcting codes 

35 such as convolut ional codes, turbo-codes, Reed-Solomon codes etc.; 
the decoding 



10 



being able to be carried out by a Viterbi algorithm, a MAP (maximum a 
posteriori) algorithm or iterated LOG-MAP etc. 

Figure 2 illustrates a general diagram of the architecture of 
the space-time decoder according to the invention. This decoder 
5 implements an iterative method making it possible to determine the 
information symbols S following a firm decision generated by a 
threshold detector 10. However, the space-time decoder according to 
the invention also generates information symbols following a soft 
decision capable of being injected into the soft input decoder 6, the 

10 injected signal being the signal obtained during the last iteration. 
The architecture of the decoder according to the invention chiefly 
involves two modules B and D injecting their output signals into a 
subtractor 9. The subtractor 9 generates an information symbol vector 
following a soft decision, these symbols being then detected by the 

15 threshold detector 10 in order to obtain symbols estimated by firm 
decision . 

The two modules B and D represent matrix products of their 
input signals by the matrices B' and D f , the index / expressing the 
iteration in progress. The module B receives at its input the signal 
20 Y from the antennae 8. The product of the matrix BL and Y is a signal 

r' from which is subtracted the signal Z from the matrix product of 
the matrix D* and the signals 1 " 1 . The matrices B' and D' are such that: 
by taking up the notations in Figure 2, we arrive at: 

7 l =B ( Y 

25 

30 Noting that in the detection of the symbol Si of the vector s, 

the i th component of z e must not be subtracted from 7 € , the following 
constraint is imposed on D': 



35 



D* = 0, Vl</<« 
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i.e. a zero diagonal in order not to subtract the useful signal 

The detector generates t 
expressed in the iteration t by 



The detector generates the firm decision S from 5 which is 



5 s l = B e • Y - D ( ' s' A = B ( • (H s+w)-D 1 ' s A 

The first step to be implemented in the iterative decoder 
consists of determining the matrices B f and D' such that the mean 
quadratic error at the threshold detector input is as small as 
10 possible. It is defined at iteration 1 by the quantity: 

MSE'(B f , D ( ) = ||s' -sf 

Minimizing the following expression means that B l and D* verify 



15 



dMSE* (B e , Z>' ) d MSE e (B e , Z>< ) 



= 0 



20 



The solving of these equations makes it possible to obtain, at 
iteration f as a function of the iteration f-1, the following 
matrices B' and D': 



i 



&s ) 



\<i<MT 



.C H V l 
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D =B<.C Pf _ i -Diag 



1 



p}-\ A f + 



-s J 



\<i<MT 



With: 
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v = 



x -CC H + Id N 



and 



A l = diag( C H V' C); 

where E s is the mean energy of the QAM constellation and N 0 is the 
noise variance. 

10 The form of the matrix D' is intuitively satisfactory. In fact, 

if S e ~ l =5 such that p e ~ l =z\, then the D' output reproduces the inter- 
symbol component for each symbol s i7 V 1 < I < M . 

In a more general manner, // 1 = 1 indicates the confidence there 
is in the quality of the estimated S*' 1 . If is not reliable then 

15 // 1 = 1 will be low and consequently a lower weighting will be applied 
to the estimator of the inter-symbol interference subtracted from 
r . On the other hand, if S iA is an excellent estimate of s, then 
Ps l — > 1 and almost all the inter-symbol interference is subtracted 
from . It should also be noted that at the first passage (=1, 
20 p* s ~ l = pQ=0 and as there is no estimated 5 M =5° still available, B' is 
the adapted filter. 

In these equations, p' designates the value, at iteration of 

the correlation between the symbols detected at iteration S, and 
the symbols actually transmitted. This correlation is therefore 

25 



<_ E <sA) 
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where E s is the mean energy of a symbol 

It can be noted that the matrices B f and D' have very different 

roles : 
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The matrix B f is a signal-formation matrix. At iteration 
zero, it is noted that B' is the linear decoder which minimizes the 
mean square error. When the correlation becomes great (tends towards 

1) , the role of the matrix B' becomes marginal. 
5 - The role of the matrix D' is to subtract from the current 

information symbol the interference due to the other information 
symbols. Its role is marginal at iteration zero, but as confidence in 

the detected symbols S* grows, its role becomes determinant. 

A means of estimating the correlation coefficient will now be 
10 described. 



t E( S J[) 

Calculation of the correlation p = requires 

E 5 

calculation of the signal to interference ratio at iteration 
SINR/. This signal to interference noise ratio is first calculated as 
15 a function of the correlation at the preceding step f-1. 
1. Calculation of SINR' 
It is shown that: 



SINR* = 



-1 



A'-J-eM') J 0-</> M > 2 ) 



20 



with 



1 - l-fc") 1 1 = N-E 5 = N 
& C ' C N 0 M-N Q 

25 

S '\ ~~J~dt , the integral exponential. 



2. Calculation of p f 



30 The calculation of p' is done in several steps: 



a - Let p° = 0 and / = 1 



b - Calculate the signal to interference (+ noise) ratio at the 
threshold detector input using the formula: 



c - Calculate the symbol error probability Pr at the threshold 
detector input from the SINR' assuming that the total noise is 
Gaussian and using the formula corresponding to the constellation, 
for example in the case of an m-PSK type modulation, the following 
formula can be used: 



This calculation results from an approximation to a strong signal to 
noise ratio for the symbol error probabilities of the m-PSK 
modulations associated with a threshold detector symbol by symbol in 
the presence of an additive white Gaussian noise type channel. This 
approximation is linked to the definition of the equivalent model. 

d - Calculate the expression of p f at the subtractor output assuming 
that when there is an error, then the threshold detector detects one 
of the closest neighbours to the symbol transmitted. In the case of 
an m-PSK, this results in: 



e - Increment f ( f = f +1) and return to step b. 

As seen in Figure 2, the architecture of the decoder according 
to the invention shows soft decisions at the output of the subtractor 
9. These decisions are exploited at the final iteration. They can be 
injected into a soft input decoder 6. 



SINK* = 




Pr 
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By calling the final iteration L, it is possible to approximate 
S 1 by the following equation: 



S = K S + noise 

5 

Where constant K > 0 is known and S is the column vector of the 
symbols with MT components. 

If the soft or weighted input decoder is a Viterbi algorithm, 

it will only have to minimize the euclidian distance between s; L and 
10 KS on all the error-correcting code words. 

Figure 3 shows a general block diagram of the decoder according 
to the invention. The decoding module 11 comprises the elements B, D, 
9 and 10 of Figure 2. In fact the vector Y at the input, and a soft 

output s 1 (output of the subtractor 9) as well as the output 
15 (output of the threshold detector 10) are found. It shows a module 12 
for calculation of the coefficients of the matrices B and D. This 
module 12 makes it possible to calculate, at each iteration or 
overall before the start of the iterations, the correlation 
coefficient and the coefficients hf and df of the matrices B and D. 
20 This module 12 receives at the input the matrix H (estimation of the 
transmission channel) and the variance N 0 . It can generate at the 
output the correlation coefficient which can be used to interrupt the 
iterations or for any other use. 

Of course the invention is not limited to the examples which 
25 have just been described and numerous of adjustments can be to these 
examples without going beyond the scope of the invention. 



